x=-1.1:0.01:1.1;
y=-0.2:0.01:1.1;
[X Y]=meshgrid(x,y);
Z=rosenbrock(X,Y);
contour(X,Y,Z,1000);
hold on;
plot (1,1,'r*')